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It has been revealed by mean-field theories and computer simulations that the nature of the 
collapse transition of a polymer is infiuenced by its bending stiffness e^. In two dimensions, a recent 
analytical work demonstrated that the collapse transition of a partially directed lattice polymer is 
always first-order as long as eb is positive [H. Zhou et al, Phys. Rev. Lett. 97, 158302 (2006)]. Here 
we employ Monte Carlo simulation to investigate systematically the effect of bending stiffness on 
the static properties of a 2D lattice polymer. The system's phase-diagram at zero force is obtained. 
. . . Depending on and the temperature T, the polymer can be in one of three phases: crystal, 

00 ' disordered globule, or swollen coil. The crystal-globule transition is discontinuous, the globule- 

, coil transition is continuous. At moderate or high values of th the intermediate globular phase 

■ disappears and the polymer has only a discontinuous crystal-coil transition. When an external 
force is applied, the force-induced collapse transition will either be continuous or discontinuous, 

(~] , depending on whether the polymer is originally in the globular or the crystal phase at zero force. 

■ The simulation results also demonstrate an interesting scaling behavior of the polymer at the force- 
induced globule-coil transition. 
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I. INTRODUCTION 



■ I ' ^ A polymer is a linear chain of monomersi. The configuration of a long polymer in solution is influenced by three 
'"^ . different types of interactions. Firstly, there are complicated interactions between monomers such as hydrogen- 
^ ■ bonding, weak van der Vaals attraction, and electrostatic interactions. Monomer-monomer interactions can also be 
O mediated by solvent molecules, e.g., the attraction between two negatively charged monomers due to the mediation 
I ^ I of a multi-valent positive charge (Mg^+ , for example) . fn theoretical models of homologous polymers, the monomer- 
monomer interactions can be represented by a potential energy U{ri,rj) between two monomers at position and 
T-H ' position Tj. Monomer-monomer interactions can be either attractive or repulsive. If the total interaction between a 
^ [ pair of monomers is attractive, it tends to bring different parts of the polymer together. At low temperatures, these 
^ . monomer-monomer attractive interactions result in the formation of compacted (globular) polymer configurations. 
' The volume of the polymer shrinks as much as possible to minimize the total monomer-monomer contacting energy. 
. In such a globular configuration, the radius of gyration Rg of the polymer scales with the polymer length N as 
i?g oc N^/'^, with d being the spatial dimension^. (For polymers in solution, d — 3; for polymers attached on the 
membrane of mobile lipids^., d — 2). 
qIq ' In the solution, solvent molecules collide frequently with the polymer due to thermal agitation. These collisions 
, tend to make the polymer to take disordered and swollen coil configurations. At high temperatures, thermal effects 
■ win over monomer-monomer attractions. Consequently the polymer will take a randomly coiled shape, whose radius 
of gyration scales with the polymer length N as Rg oc N", where the scaling exponent j/ « 3/((i+ 2) according to the 
Flory mean-field theory2i4i& and 2D exact calculations^. 
5-H Between the above-mentioned low-temperature globular phase and high-temperature swollen coil phase, there exists 

also another critical phase when the temperature is at a critical point Tg, the so-called theta-point^. At Tg, the polymer 
achieves a delicate balance among monomer-monomer attractions, excluded- volume repulsions, and configurational 
entropy. The polymer's gyration radius follow another new scaling law Rg oc N'^" , with the critical exponent i^cr — 4/7 
for 2D polymers^ii and fcr — 1/2 for 3D polymers^iii. 

The third type of interactions are external forces. In recent single- molecule manipulation experiments, external 
forces on the order of piconewtons (10~^^ N) can be applied on the ends of a polymer such as DNAS. When the 
external force field is sufficiently weak, the configuration of a self- attracting polymer is not affected by the force, since 
the attractive interaction between monomers is stronger to keep itself remain compacted. On the other hand, when 
the external force is sufficiently large, the polymer will be elongated considerably along the force direction, and the 
end-to-end distance of the polymer chain scales linearly with polymer length N. 

The overall shape of a self-attracting polymer can be dramatically changed by changing environmental temperature 
and/or external force. The globule-coil structural transition is a fundamental issue in polymer physics. The nature of 
this transition has been studied extensively in the past several decades by experimental, theoretical, and simulation 
approachea^'^ i^'^^i^^'^^i^^'^^'^^i-'^^'^^i"'^^'^^ . Certain degree of consensus has been achieved on various aspects of the poly- 
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mer globule-coil transition, but there are still some unresolved important issues. For a flexible self-attractive polymer, 
mean-field theories predicted a second-order continuous globule-coil phase-transition as the ambient temperature is 
elevated. This prediction was confirmed by more recent analytical calculation o^°i^^'^^ in 2D; however, there are still 
some controversies in the literature concerning Monte Carlo simulations in 3D (for example, Refs. 23ll2^ believed 



that the 3D globule-coil transition is a first-order phase-transition in the thermodynamic limit). The force-induced 
globule-coil transition is found to be first-order in 3D and to be second-order in 2D ^°i^^ . 

The phase behavior of a semifiexible polymer is even more complex. For a semiflexible polymer at very low 
temperature, the arrangement of its segments in a globular conformation may be highly ordered. This is driven 
by the desire of minimizing the total bending energy. The polymer are therefore in a compact crystal phase. As 
temperature increases, this crystalline order may be destroyed so as to gain configurational entropy, but the polymer 
is still in a (disordered) globular form. This crystal-globule transition is predicted to be a first-order phase-transition 
both by exactly solvable models^i^ and by 3D mean-field theory and Monte-Carlo simulation studiesi^iiii^^. When 
the temperature is further increased to the 9 temperature, the polymer globular conformation will change again to a 
swollen coil conformation. This globule-coil transition is believed to be second-orde r^i-'^^'^'^i^^'^^i^^ . The 9 temperature 
appears to be insensitive to the polymer's bending stiffnesa^^ii^'^. However, if the bending stiffness of the polymer 
is sufficiently large, the solid-globule and globule-coil transitions may occur at the same temperature, resulting in a 
single first-order solid-coil structural transitioni^. 

The force-induced collapse-transition have also been studied in recent years. For 3D polymers, the transition is 
shown to be first-order no matter whether the polymer is ficxible or semiflexibl o^^i^^ . For 2D ficxible polymer, this 
transition is believed to be second-orde r^°i^^i'^° ; while the results of an exactly solvable model^ suggest that the order 
will change to be first-order when the polymer is semifiexible. 

In this work, through a series of Monte Carlo simulations we aim at a comprehensive understanding of the phase 
behavior of a 2D semifiexible polymer under the action of temperature and external force. We pay special attention on 
the effect of bending stiffness. First we obtain the phase-diagram of the polymer at zero force using temperature T and 
bending stiffness Ch as two control parameters. At low Cb, the polymer can reside in three phases: crystal, disordered 
globule, and swollen coil. The crystal-globule transition is first-order and the globule-coil transition is second-order. 
When eb exceeds some threshold value, however, the intermediate globular phase of the polymer disappears. The 
polymer then has only one crystal-coil transition, which is first-order. After the zero- force phase-diagram of the 
polymer is known, we continue to study its stretching behavior. We find that if at zero-force the polymer is in the 
crystal phase, the applied force will cause a first-order crystal-coil transition. But if the polymer is originally in the 
disordered globule phase, the force-induced collapse transition will be continuous. In the later case, the critical scaling 
behavior of the extension of the polymer with chain length at the transition point is obtained by MC simulation. To 
predict this critical scaling behavior analytically requires further work. With these results, the present simulation 
work enriches our understanding of the structural properties of 2D self-attractive polymers. It contributes to a full 
characterization of semiflexible polymers both in 3D and on a surface. 

The paper is organized as follows: The next section introduces the model system; and the Monte Carlo simulation 
method is described in Sec. lIIII Section lTVl and Sec.|V]focus, respectively, on temperature- and force-induced structural 
transitions. We conclude this work in Sec. IVTl 



II. THE MODEL 



We study by MC simulations the phase behaviors of a 2D self-avoiding lattice polymer. The polymer is on a square 
lattice with edge length a. The bond length between two consecutive monomers i and j -I- 1 of the polymer also equals 
to a. One end of the polymer is fixed; the other end is free (without external force) or is being stretched by an external 
force / along the cc-axis of the surface (see Fig. [T]). 

If two monomers i and j of the polymer (with |j — i| > 2) occupy nearest-neighboring sites on the square lattice, 
there is an attractive energy of magnitude e. If in a configuration C there are a total number of Nc such contacts, 
then the total contacting energy is —Nee. For semifiexible polymers, there are also bending energies. In the lattice 
polymer, whenever the configuration C makes a turn (with three consecutive monomers not staying on a rectilinear 
line) there is an energy penalty eb- Denote the total number of turns in the configuration C as iVb, the total bending 
energy is then N^eh- The energy contributed by the external force is equal to —fx, where x is the a;-axis component 
of the end-to-end distance vector of the polymer. For a given configuration C, the total configurational energy is 
expressed as 



E{C) = -N,e + NbEb - fx 



(1) 
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FIG. 1: A polymer on a 2D square lattice with lattice constant o. The configuration of the polymer fluctuates over time under 
the joint action of monomer-monomer contacting interactions, bending energy, external force, and thermal energy. 



The partition function of the polymer is 



E 



exp 



(2) 



where /cg is the Boltzmann constant and T is the environmental temperature. 

In the following MC simulations, energy unit is set to be e, length unit is set to be a, fc^ is set to be unity, and force 
unit is set to be e/a. In the absence of external force, the 3D version of this lattice-polymer was previously studied 
by mean- field theory and by MC simulations in Ref. [lj|. 



III. MONTE CARLO SIMULATION METHOD 



Many different Monte Carlo procedures are documented in the literature (see, e.g., Refs. [illllllllllSiill^) to 
sample configurations of a polymer according to the Boltzmann distribution 



exp{-E{C)/{kBT)) 



Z 



(3) 



where Z is the total partition function of the polymer as given by Eg. The present MC procedure, which 

is based on the idea of importance sampling (see, e.g., textbook [37j|), is inspired mainly by the earlier work of 
Refs. [l^l3lll3^l33l | . In sampling polymer configurations, the transition probability T{fi v) for the polymer to evolve 
from an old configuration /i to a new configuration v is given by^^: 



(4) 



In this equation, — > v) is the probability of proposing a particular updated configuration v from the old configu- 
ration /i, and A{^, — > v) is the probability of accepting this updates. To ensure detailed balance, — > v) is chosen 
such that 



(5) 



where ri{fi —>■ i>) = g{v — > fi)/g{ii v). 

In this paper, we use six different types of elementary updating rules for the polymer's configuration (Fig. [2]). These 
updating rules are described in some detail here. 
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FIG. 2: Six different types of elementary configurational updates of the present MC simulation method. 



rule-1 Bonds i — 1, i, and i + 1 are parallel to each other [Fig. [21^ a)]. After the update, a hairpin is formed right after 
bond i — I, with bond i + 1 being at the head of the hairpin. The polymer segment from bond z + 3 to bond N 
reptates along its old contour and the free end of the polymer shrinks by two bonds. 

rule-2 Bond i is perpendicular to bonds i — 1 and i + 1 [Fig. (Hb)]. After the update, i becomes parallel to bond i — 1, 
while bond i + 1 is perpendicular to both bond i — 1 and bond i. 

rule-3 Bonds i ~ 1, i, and i + 1 forms a hairpin, with bond i being at the head [Fig. HIc)]. After the update, with 
probability one-half, the stem length of this hairpin increases or decreases by one bond length. This stem length 
change is achieved by reptation, with the free end of the polymer draws back or stretches out by two bonds. If 
the free end of the polymer stretches out, the orientations of the last two bonds are randomly and independently 
assigned. 

rule-4 Bond z + 1 is perpendicular to both bonds i ~ 1 and i [Fig. [2{d)] . With probability one- half, after the update 
bond i becomes perpendicular to both bonds i — 1 and j -I- 1. With the remaining probability one-half, after 
the update bond i + 1 becomes the head of a hairpin while the polymer segment from bond i + 3 to bond N 
reptates along its old contour, and the free end of the polymer shrinks by two bonds. 

rule-5 Bond i — 1 is perpendicular to both bonds i and i + 1 [Fig. [Ife)]. After the update, bond i + 1 becomes the 
head of a hairpin, while the polymer segment from bond i -I- 3 to reptates along its old contour and the free 
end of the polymer shrinks by two bonds. 

rule-6 The polymer segment from bond i -I- 1 to bond N is rotated as a whole with respect to bond i clockwise or 
counter-clockwise by an angle of tt/2 or tt [Fig. [DJf)]. This elementary update causes a large change in the 
chain's configuration. 

We have taken special care to ensure that during these MC updating processes, (1) no two monomers occupy the 

same site, (2) detailed balance is not violated. To ensure detailed balance, the ratios r]{fi v) of Eq. ([5]) are carefully 

determined and given in Fig.[2l (When the chosen bonds of Fig.[2]are located at the ends of the polymer, the — > v) 

values will be different from the values shown in the figure.) With the present MC procedure, we are able to simulate 

polymer lengths up to A^ = 3, 000, which is much longer than the lengths of several hundreds used in some previous 
MC simulationsiiiSliSS^ 
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IV. TEMPERATURE-INDUCED COLLAPSE TRANSITIONS 



We begin with the temperature-induced coUapse transitions at zero external force (/ = 0). We wiU obtain numeri- 
cally the phase diagram of the polymer in the plane spanned by temperature T and bending stiffness eb- 



A. Flexible polymers (tb = 0) 



First we focus on flexible chains. The equilibrium properties of a flexible lattice polymer are well-known, therefore 
this system can be used to check the validity of our MC simulation method. Let us denote by Rn the end-to-end 
distance of the polymer of N bonds. (In our simulations, we use Rn instead of the gyration radius to characterize 
the structural properties of a polymer. R^ behaves in the same way as the gyration radius^ and it can be computed 
with less statistical variance^^). As mentioned in Sec. HI there is the following scaling relationship between R^ and 
the polymer size N: 

{RD-N^". (6) 

where (■ ■ • ) means thermal average. At low temperatures, the polymer is in the compacted globular form, and the 
scaling exponent v approaches 1/2 in 2D^ '^^i'^° ; at high temperatures, the configuration of the polymer is reminiscent of 
a self- avoiding random walk, and v approaches the other limiting value of 3/4 [6]. At the globule-coil phase transition 
temperature Te, the polymer is in a critical state, and the scaling exponent v takes yet another value v = Vct = 4/7 
. The scaling exponent v of Eq. ^ in general may have a weak dependence on the polymer length TV— i^i^i^i^i^. 
However, at the critical state, v will be independent of N . If we plot (i?^)/iV^'^", then curves for different lengths 
N should intersect at the temperature. 

Figure [3] shows the relationship between [R^j^) /N^/"^ and temperature for chains of different lengths N. Different 
curves indeed intersect at the same point when T — Tg k, 1.50. Therefore, the present MC method is able to correctly 
predict the critical exponent of Vcr = 4/7 for a 2D flexible polymer. The predicted Tg value is also in agreement with 
the value of Tg = 1.54 ± 0.05 as given by Ref. 0. 

When there is no external force, the specific heat c{T) of the system as a function of T can be obtained through 
c(T) = 0^{{E'^) - {Ey)/N and is shown in Fig.H As the chain length N increases, surface effect (which reduces the 
total number of contacts) of the globular phase becomes less important, therefore the temperature Tp at which c(T) 
reaches maximum increases with N. The MC simulation data suggest that Tp approaches the critical temperature Tg 
according to the following formula Tp ^Tg[l- {co/Ny] with cq = 3.0 ± 1.0 and s = 0.44 ± 0.04. The peak value c{Tp) 




FIG. 3: (Color Online) The relationship between (R%) /N^^'' and temperature T for a flexible lattice polymer of different 
lengths A''. Different curves intersect at the 9 temperature of Tg ~ 1.50. 
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FIG. 4: (Color Online) The specific heat c{T) for a 2D flexible lattice polymer of different lengths A''. 
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FIG. 5: The relation between the peak values c{Tp) of the specific heat and chain length for a flexible polymer (eb ~ 0). The 
solid line is the fitting curve c(Tp) = bN" with b = 0.179 ± 0.017 and a = 0.166 ± 0.015. 



of the specific heat diverges with chain length N according to c{Tp) oc N" with a = 0.166 ± 0.015 (see Fig. [5]). From 
this finite-size scaling behavior we estimate the crossover exgonent-4''' = 1/(2 — a) to be i/i « 0.545. This value is in 
agreement with the value of (j) = 0.52 ± 0.07 given by Ref. |46| . 

Taking together, the simulation results of this subsection, in confirmation to various earlier efforts, predict a 
continuous globule-coil phase-transition for a self-attracting fiexible 2D polymer. 
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FIG. 6: Typical configurations of a semifiexible polymer in the crystal (a), disordered globule (b), and extended random coil 
(c) phase. 



B. Semifiexible polymers 



When bending stiffness eb of the polymer is positive, the polymer's configuration is also influenced by bending 
energies. At very low temperatures entropic effect is not important, then the polymer will favor those compact 
configurations which have the maximal number of monomer- monomer contacts and the minimal number of bends. 
These are highly ordered crystal configurations as illustrated in Fig. [6ja). fn such a crystal configuration^- ^^'"'^ long 
segments of the polymer stack onto each other, and the total number of bends iVb scales only sub-linearly with chain 
length N. Let us define two order parameters nb (the density of bends) and ric (the density of contacts) as 



rib 



N 
N 



(7) 

(8) 



For a polymer in a crystal conformation, we have limjv^oo nb = and limjv_»oo nc of order unity. The overall shape 
of a polymer crystal may be a rectangle rather than a square, the aspect ratio of which depending on eb- When 




FIG. 7: (Color Online) Evolution of the configuration of a semifiexible polymer with bending stiffness eb = 0.3 and length 
N — 2, 000. The vertical axis is the number of bends A'b in the configuration, and the horizontal axis is the evolution time in 
terms of MC steps. (At time zero of the figure the system has already reached equilibrium.) The temperature is T = 0.33 (red 
stars), 0.38 (black open squares), and 0.39 (green filled squares). 
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FIG. 8: (Color Online) The relationship between {R%) /N^^'' and temperature T for a semiflexible polymer of different lengths 
A^. The bending stiffness is eb ~ 0.6. Different curves intersect at the 9 temperature of Te « f .55. 

temperature is elevated, the entropy effect becomes more and more important. At certain point the crystal order of 
the polymer will be destroyed. For Cb being small, when crystal order disappears the monomer-monomer contacting 
interactions may still be very significant. Then the polymer will take disordered globular shapes (Fig. EJ^b)). In a 
disordered globular configuration, both the order parameters Uh and ric are of order unity. When temperature is 
further increased, the globular configuration will again become unstable, and the polymer will transit to the phase of 
swollen random coil (Fig. (S^c)). 

The crystal- globule structural transition is a first-order phase-transition in the thermodynamic limit of iV ^ oo, 
similar to the solid-liquid transition of water at 0°C. This is demonstrated in Fig. [7] for the evolution of the total 
bending number A^b of a polymer with length N = 2, 000 and bending stiffness eb — 0.3. At temperature T = 0.33, A'b 
fluctuates around Nj^^"^ w 410. At T = 0.39, iVb fluctuates around the value of N^^'' « 510. While at the intermediate 

temperature T — 0.38, TVb jumps between these two values Nj^^"^ and indicating that the system is bistable. On 

the other hand, at each of these temperatures, the histogram of the number of contacts Nc has only one peak. The 
crystal-globule transition temperature increases with the bending stiffness eb of the polymer. 

When the temperature is further increased to T « 1.5, a second structural transition takes place. This globule-coil 
transition is second-order, as can be inferred from Fig. [8l This figure shows the value of (i?^)/-/V®/^ as a function 
of temperature for a semiflexible polymer of bending stiffness eb = 0.6 and different lengths N (similar results were 
obtained for eb = 0.3 and eb = 1-0). Figure [5] is qualitatively the same as Fig. [3] of a flexible polymer. The only 
quantitative difference is that the transition 9 temperature is slightly increased. At the globule-coil transition point, 
the scaling exponent of the semiflexible polymer is the same as that of a flexible polymer. Therefore, bending 
stiffness does not influence the scaling behavior of a polymer in its critical 6 state. The globule-coil transition 
temperature is also not very sensitive to the bending stiffness. 

In the crystal or the disordered globular phase, the size of the lattice polymer scales with chain length in a square- 
root law, i.e., Rj^ = ciN with ci w 0.6. The value of the prefactor ci has no obvious change at the crystal-globule 
transition temperature; it is also insensitive to the bending stiffness eb- When the temperature approaches or goes 
beyond the globule-coil transition temperature Tg, however, ci increases quickly with T (because the square- root 
scaling law no-longer holds). 

When the bending stiffness of the semiflexible polymer exceeds a threshold value of eb ~ 1-5, the intermediate 
disordered globule phase disappears. The polymer will exists only in two phases, either the low-temperature crystal 
phase or the high-temperature swollen coil phase, and the structural transition between these two phases is first- 
order. For a semiflexible polymer with eb = 5 and length N, Fig. [9] shows the relationship between the density 
of contacts ric and temperature T. As N becomes sufficiently large (e.g., N — 1,600) we notice a dramatic drop 
of ric at T « 3.0, indicating the existence of a globule-coil transition. For T close to this transition temperature, 
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FIG. 9: (Color Online) The order parameter ric as a function of temperature T for a semiflexible chain with eb ~ 5.0 and 
different lengths A*'. 




FIG. 10: (Color Online) Distribution of contact number A'^c (a) and bending number A'b (b) for a semiflexible polymer with 
bending stiffness eb ~ 5.0 and length A*' = 1,600. Different curves correspond to slightly different temperatures. 

Fig. [TUT a) and Fig. [TUT b) report the distributions of the polymer's contacting number Nc and bending number N^. 
When the temperature T is close to T = 3.0, there are two peaks in both probability distributions, indicating the 
system is bistable (both in terms of A^c and iVb) in this temperature range. The positions of these two peaks in the 
distribution of Nc and that of A^b remain almost fixed when the temperature changes, but their weights changes. 
As temperature increases, the system has higher probability to be in the coil phase. We can define the crystal-coil 
transition temperature as the temperature at which the system has equal probability to be in the crystal and the coil 
phase. This transition is between the crystal phase and the coil phase, because both the order parameters Uc and rib 
will jump in the thermodynamic limit of iV — > cx) rather than just one of them. 

We have performed the same analysis for other values of the bending stiffness eb. At zero external force, the complete 
phase-diagram of a 2D lattice polymer is shown in Fig. 1111 When eb < Cb ~ ^-^i polymer has two structural 
phase-transitions as temperature changes, a first-order crystal-globule transition followed by a higher-temperature 
continuous globule-coil transition. When e > e^, the polymer will transit directly from the crystal phase to the coil 
phase. In Ref. [l^j, it was predicted based on a partially directed 2D lattice model that, the collapse transition is 
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FIG. 11: The zero-force phase-diagram of the 2D lattice polymer. Square symbols are simulation results. The crystal-globule 
transition and the crystal-coil transition is first-order (as indicated by dashed boundary lines) in the thermodynamic limit, 
while the globule-coil transition is continuous (solid boundary line). 



first-order as long as the bending stiffness is positive. This apparent contradiction between Ref. (221 and the present 
work is easy to understand. In the partially directed polymer model, the compact phase, since it is highly ordered, 
is not the disordered globule phase but rather the crystal phase. Because of the additional constraint of partial 
directness, the entropy of the compact phase is suppressed compared to that of the disordered globule phase. The 
simulation results of Fig. [TT] also suggests that the crystal-coil transition temperature increases almost linearly with 
the bending stiffness Cbi in consistence with the earlier work of Doniach and co-workers on 3D lattice polymers^^. 



V. FORCE-INDUCED PHASE-TRANSITIONS 



With the zero-force phase-diagram Fig. [TT] being obtained for a 2D lattice polymer, we now proceed to study its 
force-induced structural transitions. Initially the polymer is in a low temperature crystal or disordered globule phase. 
An external force along the x-direction is applied to the free end of the polymer, and the x-component of the polymer's 
extension, Xjv, the contacting number Nc, the bending number N^, are recorded during the stretching. 



A. The globule-coil transition 

For a 2D lattice polymer with bending stiffness eb < 1-5, we can choose an appropriate temperature value (say 
T — 1.2) such that, at zero external force, the polymer is in the disordered globular phase. At this temperature, when 
the external force / becomes sufficiently large, the globular shape of the polymer will be destroyed. The polymer 
will be in a highly extended coil configuration. As has been predicted by several earlier theoretical wor k^^i^^ , in the 
thermodynamic limit of A'^ — > oo, a globule-coil transition will take place at certain critical force /cr(T). In agreement 
with Ref. ^25| , we find this 2D force-induced globule-coil transition is a continuous transition. 

The continuous nature of the force-induced globule-coil transition is checked by investigating the fluctuation of the 
elongation X^. Figure [T^ shows how the extension susceptibility Xf of a flexible polymer system (eb = 0) at T = 1.2 
changes with force /. The extension susceptibility is defined by 

Xj^^^^P[{Xj,)-{X^r]/N. (9) 
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FIG. 12: (Color Online) The extension susceptibility xf ^ function of external force / at T = 1.2 for a polymer with 
eh — and chain length TV. Insert shows the distributions of normalized mean extension {Xm)/N for the polymer system with 
= 1, 600 near the globule-coil transition. 




FIG. 13: (Color Online) The value of {Xn^/N^'''* as a function of external force / for a flexible polymer (eb = 0) of different 
lengths A'^. Different curves intersect at / ~ 0.295. 



Figure fT2l demonstrated that, as the chain length N increases, the peak of Xf becomes more pronounced (eventually 
it will approach infinity). This behavior is very similar to the divergence of the specific heat c{T) at the temperature- 
induced globule-coil transition (see Sec. IIV|) . By finite-size scaling and extrapolating to the case of = oo we obtain 
a transition force fcr ~ 0.3 for the system (an more precise estimate of the transition is reported below). The inset of 
Fig. [T2l shows the distribution of the extension of the system at different external forces. When the external force is 
close to the transition force fcr, the extension distribution has only one peak which moves with force. This behavior 
further confirms the continuous nature of the structural transition. 
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FIG. 14: The critical scaling exponent /icr at the globule-coil transition of Eq. pup as a function of the bending stiffness e;, at 
temperature T = 1.2. 



Under an external force, we can write down the following effective scaling relation between the mean extension of 
the polymer and the chain length N: 

(Xn)^ ^ m (10) 

with effective exponent fi. In the globule phase, the polymer is in an elliptic shape and /i in Ec^. (llOp is less than 2 but 
larger than 1. fi may increase with force / and reach a limiting value /icr at the globule-coil transition point. When / 
further increases and the polymer is in the coil phase, the extension of the polymer scales linearly with chain length 
iV, i.e., fi = 2. As in Sec. llVi we expect that, at the globule-coil transition point, the critical exponent /icr will be 
independent of chain length N. Then if {Xn)"^ /N^" is plotted as a function of /, curves corresponding to different 
chain length TV will intersect at the same point /„.. Figure [T51 demonstrates that this is indeed the case, and from this 
plot we estimate that fcr = 0.295 and //cr — 1-74 for a flexible polymer at T = 1.2. This value of fcr is in agreement 



with Ref. [25|, which predicted fcr ~ 0.30. To further validate this method of determining the critical force, we have 
used it to determine the critical force for a flexible 2D partially directed polymer at T = 0.59 as studied in Ref. [23 |. 
We obtained that the critical force /„■ = 0.45 (very close to the exactly known value 0.46) and /^cr = 1-2. 

For T = 1.2 we have determined the critical exponent for 2D lattice polymers with different bending stiffness eb- 
The critical exponent /icr is not a constant. It drops from the value of /icr — 1-74 for eb = to /icr ~ 1-55 for eb — 1.0 
(see Fig. Intuitively this observation is natural to understand. At the force- induced globule-coil phase-transition 
point, the polymer takes the shape of an elongated ellipse^^. The bending stiffness eb affects the shape of this ellipse 
and therefore affects the critical exponent /icr- On the other hand, to derive an explicit expression for the relationship 
between /z and eb is a hard task. The critical exponent /icr may also depend on temperature. 

The globule-coil transition force fcAT) decreases with temperature almost linearly and vanishes at the 0-temperature 
Tg, i.e., /cr(r) ^ {Te - T) for T < Tg (see Fig. [13. This linear relationship holds both for flexible and semiflexible 
chains. 



B. The crystal-coil transition 



At low temperature and zero force, a considerably stiff semiflexible polymer is in the crystal phase. We now 
investigate the force-induced crystal-coil transition by MC simulation. Since at / = 0, the temperature-induced 
crystal-coil transition is already a first-order transition, we expect the force-induced transition should also be first- 
order. This is confirmed by MC simulation results. Figure [THT a) shows the evolution of the total contacting number 
Nc of a polymer of length = 1, 600 and bending stiffness eb = 5.0 at T = 2.5. At the three different force values, 
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Temperature 



FIG. 15: The globule-coil transition force fci{T) as a function of temperature T. 




Evolution Time N IN 



FIG. 16: (Color Online) The time evolution (a) and probability histogram (b) of the contacting number A'c of a polymer with 
A'^ = 1, 600 and eb = 5.0. Different data sets correspond to different external force /, while the temperature is fixed at T = 2.5. 

the contacting number jumps between the coiled shapes of Nc « 300 and the crystal shapes of Nc ~ 1,000. The 
histogram of the contacting number as shown in Fig [TBTb) has two well-defined peaks for force / ~ 0.74. When 
/ deviates remarkably from this value, however, one of the peaks in the contacting number histogram disappears, 
indicating the system is either in the crystal phase or in the coil phase. 

Although the crystal-coil transition is a first-order phase-transition, if the chain length N is too small, the jump in 
the order parameter ric or rib can not be observed. This is because the correlation length at the crystal-coil transition 
may exceed the size of a small system. Figure [TTl demonstrates that, for a polymer with bending stiffness eb = 5.0 
at temperature T = 2.5, only chains with length N > 1,000 will show clear signature of a discontinuous structural 
transition. When the chain length N is equal to 400 or 500, at the crystal-coil transition force (which corresponds to 
the maximum of the extension susceptibility x/), the two peaks of the distribution of ric have the trend of merging 
into one peak. For these later small systems, the crystal-coil transition resembles a continuous phase-transition. 

According to the phase-diagram Fig. [Til at fixed temperature T, there exists a threshold bending stiffness value 
e^(T) at which the force- induced collapse transition changes from being second-order to being first-order. To determine 
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Relative Extension 



FIG. 17: (Color Online) Distribution of normalized mean extension {Xn)/N for a polymer with eb ~ 5.0 at T = 2.5. The 
external force is set to the value at which the susceptibility Xf reaches maximum, with / — 0.68, 0.67, 0.70, 0.72, 0.73 and 0.74 
for N = 400, 500, 700, 1, 000, 1, 400 and 1, 600 respectively. 

the value of precisely, however, is not easy. When the bending stiffness of the polymer is just slightly above ej, 
although an infinite chain will have a jump in the order parameters ric and rib at the collapse transition point, for an 
finite chain such jumps in Uc and Ub will be smeared out by large fluctuations of A^c and A'b at the transition point. 
For a polymer of length N less than or comparable to the correlation length of the system at the first-order collapse 
transition point, the elongation of the polymer may still follow the scaling rule Eq. (flO)) . i.e., resembles a continuous 
structural transition. At T — 1.2, the critical exponent /icr (defined in the previous subsection) as a function of Cb is 
shown in Fig. [T3]for polymers of length up to iV = 1, 000. We find that /icr first decreases with eb; but for et, > 1.0 
it increases with eb. We interpret this change of trend as follows: At a continuous collapse transition, the stiffer the 
chain is, the more compact it prefers and therefore the smaller the critical exponent /i is; on the other hand, at a 
discontinuous collapse transition, the stiffer the chain is, the more extended the average elongation of the polymer is 
and the larger the value of /Zcr is. Based on this argument, we estimate the value of ej^ as the point where the scaling 
exponent /icr reaches the minimal value. 

VI. SUMMARY AND DISCUSSION 

The qualitative properties of temperature- and force-induced collapse transitions of a 2D self-attractive semiflexible 
lattice polymer were investigated by Monte Carlo simulations. The system has three possible phases: the crystal phase, 
the disordered globular phase, and the random coil phase. The crystal- globule and the crystal-coil transitions are 
both first-order phase-transitions in the thermodynamic limit; while the globule-coil transition is a second-order phase 
transition. The disordered globular phase is absent for polymers of considerable stiffness. These simulation results 
are consistent with earlier theoreticali^ and simulatioi*i^iS^ studies on temperature-induced collapse transitions of 3D 
lattice polymers. They also confirm and extend the earlier simulation^ and analytica l^^i'^° studies on force-induced 
collapse transitions of 2D lattice polymers. 

The simulation results also suggested that, at the continuous force- induced globule-coil transition, the extension 
of the polymer scales with polymer length according to {X^) = aoN~i^ , with both oq and /icr being independent of 
chain length. Such a critical behavior is similar to the well-known critical scaling Rg oc aiN'^"' at the temperature- 
induced globule-coil transition [see Eq. ([6])], but it is yet to be confirmed by analytical calculations. According to 
the simulation results, the bending stiffness does not affect the critical exponent i>cr- However, the scaling exponent 
/icr for the force- induced globule-coil transition changes with the bending stiffness eb (see Fig. [T4|. We suggest such 
a qualitative difference should be amenable to experimental verifications^!^. A possible reason for this qualitative 
difference were also given in this paper. 
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For 2D polymers, we have shown that the temperature- and force-induced cohapse transitions have the same order. 
This may not be the case for 3D polymers. For 3D flexible polymers, the force-induced globule-coil transition is a 
first-order phase-transition^^; the order of the temperature-induced collapse transition, on the other hand, is usually 
regarded as second-order in many references (but Binder and co-authors^"^ argued that it is actually also a first-order 
phase-transition) . 
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